CMB EB and TB cross-spectrum estimation via pseudo-spectrum techniques 
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We discuss methods for estimating EB and TB spectra of the Cosmic Microwave Background 
anisotropy maps covering Umited sky area. Such odd-parity correlations are expected to vanish 
whenever parity is not broken. As this is indeed the case in the standard cosmologies, any evidence 
to the contrary would have a profound impact on our theories of the early Universe. Such correlations 
could also become a sensitive diagnostic of some particularly insidious instrumental systematics. In 
this work we introduce three different unbiased estimators based on the so-called standard and pure 
pseudo-spectrum techniques and later assess their performance by means of extensive Monte Carlo 
simulations performed for different experimental configurations. We find that a hybrid approach 
combining a pure estimate of B-mode multipoles with a standard one for ^-mode (or T) multipoles, 
leads to the smallest error bars for both EB (or TB respectively) spectra as well as for the three 
other polarization-related angular power spectra {i.e. EE, BB and TE). However, if both E and B 
multipoles are estimated using the pure technique the loss of precision for the EB spectrum is not 
larger than ^ 30%. Moreover, for the experimental configurations considered here, the statistical 
uncertainties - due to sampling variance and instrumental noise - of the pseudo-spectrum estimates 
is at most a factor ^1.4 for TT, EE and TE spectra and a factor ^ 2 for BB, TB and EB spectra, 
higher than the most optimistic Fisher estimate of the variance. 

PACS numbers: 98.80.-k; 98.70.Vc; 07.05.Kf 



I. INTRODUCTION 

A reliable and complete characterization of the polar- 
ized Cosmic Microwave Background (CMB) anisotropics 
is one of the main challenges in observational cosmology, 
targeted by a large set of on-going [T] |2] , being-deployed 
[3j-|6| or planned [71 [8] experiments. CMB (linear) polar- 
ization is completely described by two Stokes parameters, 
Q and [/, which are mapped over the celestial sphere by 
CMB experiments. In the harmonic domain, the polar- 
ization field can be described either with help of spin-2 
and spin-(— 2) or gradient, and curl, 5, multipoles. 
From the physics point of view, the gradient /curl de- 
composition is more natural as it is directly linked to the 
cosmological perturbations produced in the primordial 
universe. Indeed, for symmetry reasons, scalar pertur- 
bations can produce solely £^-mode-like polarization pat- 
terns and therefore the B component of the polarization 
field can be viewed as a direct tracer of the primordial 
gravity waves [9 and thus as a window onto the pri- 



* |julien.grain@ias.u-psud.fr| 

^ tristram@lal.in2p3.fr 

^ radek@apc.univ-paris-diderot.fr 



mordial Universe. Due to the presence of a secondary B- 
mode contribution generated by the gravitational lensing 
of the CMB photons by the Universe's large scale struc- 
ture [10^, such a picture is not fully correct, however it is 
hoped that the latter signal can be accurately subtracted 
permitting a recovery of the primordial 5-modes. 

In standard cosmology, the physics governing photon 
propagation is parity invariant resulting in vanishing odd 
parity correlations, {TB) and {EB). Nevertheless in 
the context of the next generation CMB B-mode exper- 
iments, the EB and TB spectra will play an important 
role, both for the data characterization and their scien- 
tific interpretation. This is because on the one hand, 
these odd-parity cross- spectra are comprehensive, end- 
to-end, null tests of the presence of instrumental and/or 
astrophysical systematic effects still present in the data 
(see e.g. [TT| |12]). On the other hand, as some non- 
standard cosmological mechanisms could produce non- 
vanishing odd-parity cross-spectra, their detection could 
become a smoking gun of such effects with potentially 
far-reaching consequences for our understanding of the 
Universe. Examples of such mechanisms include a pri- 
mordial stochastic magnetic field, which generates TB 
and EB correlations, if this magnetic field possesses a he- 
lical component pISHIS] . Similar effects can be obtained 
due to a rotation of the plane of linear polarization of 
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the CMB photons travehng from the last scattering sur- 
face to our detectors. This could result from either the 
Faraday rotation induced by interaction with background 
magnetic fields \TT-'20^ or interactions with pseudoscalar 
fields [21 . In all these circumstances a robust and reli- 
able EB and TB spectrum estimator could be needed to 
constraint, and potentially correct for, such effects. 

The goal of this paper is to investigate some of possi- 
ble extensions of the existing power spectrum estimation 
methods to incorporate the odd-parity power spectra. 
Hereafter we work within a framework of the pseudo- 
spectrum techniques. In this context the polarized CMB 
multipoles can be calculated with help of either the stan- 
dard or pure estimators, with the latter specifically de- 
signed to suppress (or to nearly suppress) E-io-B and 
B-to-E leakages. In this work we use a specific prescrip- 
tion for calculating those as introduced in [23 and later 
elaborated on in [2T^!25], for some alternatives see, e.g., 
^29.-.31j . The relative merit of using either the standard or 
the pure estimates of the pseudo-multipoles has been dis- 
cussed in the context of the B-mode power spectra show- 
ing that the pure approach leads to smaller error bars 
[23-25 . This can be understood as follows. The pure 
estimator by construction removes some spatial modes, 
referred to hereafter as ambiguous modes, which are the 
source of the leakage between the two polarization com- 
ponents. This would unavoidably lead to some informa- 
tion loss, refiected in an increase of the final variance 
of the estimated quantity, were not for a gain incurred 
as a result of removal the £^-mode power leaked to the 
5-spectrum and residing in the ambiguous modes. The 
latter effect offsets the former one at least as long as 
E ^ B diS indeed is expected for the CMB signal. Conse- 
quently, a significant gain in precision is usually observed 
as compared with the standard estimator, which retains 
all the modes together with the leaked power. By the 
same token, for estimating the E-mode power spectrum 
the standard approach is expected to be more efficient 
as this time extra variance due to the power of the B 
modes leaked into E is by far subdominant to that re- 
sulting from the removal of the ambiguous modes. 

The question, whether to include or not the ambigu- 
ous modes, and whether to do it consistently or only 
partially, is less straightforward to answer in the context 
of the two odd-parity power spectra estimation and has 
not been addressed to date on either a qualitative or 
quantitative level. This is this problem, which is at the 
focus of this work. To address it hereafter we construct 
and investigate performance of three alternative EB 
pseudo-spectrum estimator^ 



^ In principle, a fourth alternative is possible: using standard es- 
timation for the S-modes and the pure estimation for the E- 
modes. However, such an alternative is a priori disfavored as a 
significant amount of information on the £'-modes is lost and the 
5-modes still suffer from E-to-B leakage and therefore will not 
be considered in this work. 



• a standard approach [27l [28j involving the standard 
pseudo-harmonic coefficients for both the signal maps 
and no explicit correction for the E-to-B leakage; 

• a pure approach^ when both sets of the multipoles are 
pure; [23-25 ; and 

• a hybrid approach^ where the E multipoles are 
obtained from the standard estimator and the B signal 
from the pure estimator [23^ 2A^ . 

Temperature multipoles are always estimated using the 
standard pseudo-multipole approach, it may therefore 
appear that there are only two possible estimators for 
the TB spectra corresponding to B multipoles being ei- 
ther standard or pure. In fact, the situation is some- 
what more complex due to the fact that, as discussed in 
an unbiased estimation of the TB spectrum 
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requires treating it together with the corresponding TE 
spectrum. Consequently, we arrive again at three types 
of estimators of possible interest, which mirror the case 
of the EB estimators, and are therefore defined by the 
type of multipoles used for E and B signals. We will 
refer to them as standard^ pure^ and hybrid estimators. 

Our paper is organized as follows. In Section [TTj we 
present our conventions and notation for describing po- 



larization fields on the sphere. Section |III] is devoted to 
the defini tion of the different pseudo- Q estim ators of the 
EB (Sec. |IIIB| ) & TB spectrum (Sec. |IIIC[ ), with a par- 
ticular emphasis on the computation of the mode-mode 
coupling matrices. The inffuence of the sky apodization 
and pixelization scheme is discussed in Section [IV| Fi- 



nally, the performances of each estimator is discussed in 
Section [V| by quantifying the uncertainties of the recon- 
structed EB and TB spectra considered as null tests. In 
that Section, we assume a typical sky-coverage and noise 
of suborbital experiments, as inspired by the on- going 
EBEX experiment [32], and consider a series of mock ob- 
servations with progressively more realistic features. A 
simple Fisher analysis estimate of the power spectra un- 
certainties, which provides a lower limit of the variance of 
the estimators, is used as a theoretical benchmark. Our 
conclusions are presented in the last Section (Section [Vl| 
of this paper. 

Most of the technical details are provided in the ap- 
pendices of this paper. Appendix |A] summarizes the 
main properties of Gaunt integrals and Wigner-3j sym- 
bols used in our derivation of the different mode-mode 
coupling matrices, while the details of their computation 
and explicit formulae are provided in Appendix [B] In Ap- 
pendix [C| we provide the noise contribution expected to 
bias the pseudo-C^'s for the three formalisms. 

In the following of this paper, the term mask will sys- 
tematically denote the binary, pixelized maps assuming 
two values, or 1, corresponding to the observed (or kept 
in the analysis) or unobserved pixels. The pixel weights 
which can assume any non-zero value will be referred to 
sky apodization or window function. 
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II. CONVENTION AND NOTATION 

A. Spin fields on the sphere 

Any spin-5 field on the sphere satisfying = 
(in the following, s is always assumed to be > 0) can be 
decomposed using an E and B basis, 



im 



(1) 



Im 



where igl^rn stands for the spin- weighted spherical har- 
monics obtained by applying the spin-raising and spin- 
lowering operators {d and B) on the standard (spin-0) 
spherical harmonics. 

sYim = ^d'Y^^ and sY^m = B'Y^^, (2) 

with 



{i-s)V 



We remind that sY^^ = (— l)^+^_sY£(_^). For the spe- 
cial case of a spin-0 and real valued field (like the CMB 
temperature field), the B component is vanishing. More- 
over, this also applies to any spin- 5 field built from a 
spin-0, real- valued field by use of the spin-raising opera- 
tor (i.e. = d^^o)^ as can be directly seen by applying 
the spin-raising or lowering operators on the harmonic 
representation of the field as in Eq. ([T]) and using Eq. 

CMB anisotropics are described by three fields: a spin- 
field, corresponding to temperature anisotropics T(n), 
and spin- (±2) fields describing the polarization denoted 
±2P- In terms of the Stokes parameters, the two polar- 
ization fields are given by ±2P = (Q^iU). The multipole 
decompositions of such anisotropics are given by, 



(3) 



±2P(n) = - V (af, ± iafj ±2l^£m(n). (4) 



im 



Conversely, at least as long as CMB is mapped over the 
full sky, the T, E and B multipoles can be reconstructed 
from the Stokes parameters by applying the following 
scalar product. 



/ T{n)Y;^{n)dn, 
c^fm^-\l [2P2Y;^^-2P-2Y;jdn, (5) 

^fm f [2P2YL --2 P-2YIJ dfi. 

Hereafter we will use a^, (X = T, E^ B) to denote, the 
full-sky multipoles as in Eq. ([5| and referred to them as 
the CMB multipoles and use to denote partial-sky, 
pseudo-multipoles defined in the next Section. 



III. PSEUDO-SPECTRUM ESTIMATORS FOR 
THE EB & TB CROSS-CORRELATION 

A. Standard and pure pseudo-multipoles 

1. Temperature case 

For temperature anisotropics, the pseudo-multipoles 
are defined as 



An, 



(6) 



where W^^^ is a window function applied to the temper- 
ature map. Such pseudo-multipoles are given by a con- 
volution of the CMB temperature multipoles aj^ with 
multipoles of the window function. 



J-im — J^im;i'm'^i'm'' 



(7) 



An exact expression for the convolution kernel K^^.i/^i 
has b een derived in [271 [Ml ESI |34| and is recalled in Eq. 

(|bT4|). 



2. Polarized case: standard pseudo-multipoles 

In the case of the polarization, pseudo-multipoles can 
be introduced most straightforwardly by directly adapt- 
ing the definition from the temperature case, Eq. ([6|, 



{std) 



B 



{std) _ ^ 



\l W[2P2Yl^^ -2P-2YlJdn, (8) 
j W [aP 2^/^ - -2P -2YIJ dn. (9) 



We will refer to an approach based on these definitions as 
the standard estimator. The resulting standard pseudo- 
multipoles are then related to the CMB polarization mul- 
tipoles by. 



^im — y^im.i'm'^i'm' ^ im,tm'^i'r 



B 



(std) 
im 



_ fr^diag B 

— [^im.i'm'^i'n 



^-^iml'm'^i'm' 



, (10) 

• (11) 



The convolution kernels, K^^^^ and K^^^ ^ are functions 
of the Gaunt integrals and their explicit expression are 
given in Eq. (B7) of App. B 1 



3. Polarized case: pure pseudo-multipoles 

As emphasized in [35], the pseudo-multipoles as de- 
fined above suffer from E-to-B leakage. Following [35] . 
the polarization field on a partial sky can be split into 
three subspaces: the pure £^-modes, the pure 5-modes 
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and the ambiguous modes, which contain simultaneously 
information about E and B modes. These are these 
last modes, which are the source of the E-io-B leakage 
due to impartial sky coverage. Their presence is due to 
the fact that on the cut-sky the spin harmonics, ±2y£mi 
are not anymore orthogonal and any alternative basis 
constructed out of their linear combinations will include 
some which will be neither of the B- or E- type. 

Because, for the CMB the power of the E signal is 
much bigger than that of 5, this leakage is a domi- 
nant factor for the 5-modes estimation giving rise to a 
so called E-to-B leakage problem. Nevertheless, by the 
same token a certain amount of B-mode power is also 
present in the E pseudo-multipoles, which occasionally 
may also become problematic. 

The E-to-B leakage problem can be resolved if the 
ambiguous modes are identified and excluded from the 
spectrum estimation procedure or by using some (incom- 
plete) set of basis functions, which by construction are 
orthogonal to such modes. Following [23] such a basis 
can be constructed with help of the spherical harmon- 
ics and appropriate window functions W and Wi = dW 
both vanishing at the boundaries of the observed sky. 
Denoting these basis functions as. 



(12) 
(13) 



we can write now expressions for pseudo-multipoles, 
which are free of any E/B leakage as. 



^{pure) 



B 



{pure) 



'2L 
2/. 



P V 



p viW)-" 

2^ -2yfm. 



2P2yim 



P-9Y, 



df{;i4) 

df{15) 



As in the standard case these multipoles can be related 
to the actual CMB multipoles as. 



E, 



{pure) 
im 



Tj-diag E ^AZjoff B 



B 



{pure) _ \rrdiag B _^zroff E 

im ~ Y^im.i'm'^t'm' ^ tm ,t' m' ^ t' m' 



,(16) 
.(17) 



We note that in principle by construction i^^^^-^) should 
be identically zero. This would be however only the case 
were the cut-sky a sole source of the leakage. In prac- 
tice other sources of residual leakage exist, most notably 
due to sky pixelization effect^ and it is prudent and use- 
ful to correct for such effects at least on average, what 
can be achieved with help of an appropriately calculated 



^ More particularly, the Dirichlet and Neumann boundary condi- 
tions can never be completely fulfilled by W because of pixeliza- 
tion. 



off-diagonal block. Such off-diagonal blocks, though not 
precisely zero, have been shown to be typically much 
smaller than the diagonal blocks for the pure estimates. 



I.e.. 



B^ 



off 



diag 
im,tm' 



and, more importantly, than 



the off-diagonal blocks of the standard approach, i.e., 

HeZ'm' « K^Z'm' ESII^- This is this latter condi- 
tion, which ensures that the pure estimator typically re- 
sults in significant improvements over the standard one at 
least as far as B-mode power spectra are concerned. This 
emphasizes the fact that an estimator does not have to be 
strictly pure, i.e., H^^^ / 0, to provide good estimates of 
the B modes. This indeed turns out to be the case with 
so called optimized apodizations [24 , which lead to the 
smallest overall uncertainty of the estimated BB spec- 
tra even if not strictly fulfilling the required boundary 
conditions. 

The pure pseudo-multipoles, Eqs. (14) 



15L can be 



related to the standard pseudo-multipoles as follows. 



{pure) 
im 



E2,im + '^ITT^El/m + ^0,£m, (18) 



^(pnre) ^ B2/m + 2— ^^i^^^ + — — Bq^^^. (19) 

Here, Es/m{Bs/m)^ s = 0,1,2 is a spin-s multipole of 
type E (b) of' a field, W^_, 2P , where Ws = d' W 
are derivatives of the scalar window, W. It is easy to 



show that E2/ 



_ T^{std) 



im 



and B2/m B 



{std) 



while the 



two remaining terms, called the counter-terms, are there 
to cancel the ambiguous modes present in the standard 
term, what they do prefectly whenever the window ful- 
fills the required conditions. Each of the spin-5 pseudo- 
multipoles is related to the CMB multipoles via spin- 
dependent, mode-mode coupling matrices as in Eqs. (10) 



and (11). Combined together, those lead to expressions 
for the convolution kernels, H^'^^^ and H^^^ ^ introduced 
in Eqs. (16) & (17). The explicit forms for all those ma- 



trices can be found in App. B 1, Eqs. (B7) and (B7) 



4. Coupling between aim and pseudo-aim 

It is instructive to briefly comment on how the pseudo- 
multipoles are coupled to the CMB and the window func- 
tions multipoles. The mode-mode coupling corresponds 
to the coupling of three angular momenta: the momen- 
tum (-^, m) of either the E or B pseudo-multipoles, the 
momentum {^' ^ m') of either the E oi B CMB multipoles 
and the momentum {I" ^ m") of either the E 01 B compo- 
nents of each spin-5 window functions. The total number 
of such angular momentum couplings is eight as any an- 
gular momentum can be of either E- or 5-type. Each of 
these couplings can be classified according to its parity: 
even or odd. For example, the coupling between Eirni 
aff^, and w^f^J,^,, does not vanish only for even values 
of (£ + £'+ £''), as it is made of the quantities (see 
App. [A|. We call such a coupling even. In contrast, the 

coupling between Eirn^ ^f'm' ^^f^i"m"'> ^^^^ of "^^e 
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G~ quantities, is an odd coupling as it is non-zero only 
for odd values of {i ^ i' ^ i"). From this perspective, 
the pseudo-multipoles can be written using the following 
matrix notation 



odd 



_AnOdd 
r,even 



X 



W 



W 



(20) 



Bem = XI X! X! ( 



' s=0 



odd 



-w. 

^even 



even 
s 

odd 



W 



(21) 



where both gf^^^ and g^/"^ depend on (^, m), {I' ,m') and 
{e',m"), and ^/f x gf^ = if they both have the 
same arguments. This classification is helpful during the 
computation of the pseudo-power spectra as any term 
involving the product of an odd coupling times an even 
coupling will vanish. The parity nature of each coupling 
is summarized in Table [H 





E (E) 


E (B) 


B (E) 


B (B) 


Eirn 


even 


odd 


odd 


even 




odd 


even 


even 


odd 



TABLE I. Parity classification of the coupling between 
pseudo-multipoles, CMB multipoles and window function 
multipoles. 



where (x) and (y) stands for the standard and the pure 
pseudo-multipoles respectively. Because of the limited 
sky-coverage and pixel dependent weights, these estima- 
tors are biased as their average over CMB realizations in- 
volves a mixing between different £-modes and between 
different polarization type. Those two classes of mix- 
ing are encoded in the mode-mode coupling matrices 
M^^f ^ unbiased estimator, C["^"\ is com- 

puted by solving the following linear system 



EE, EE 
BB,EE 



j^EE,BB 
ilfBB,BB 
j^EB,BB 



BB,EB 



EE,EB \ / ^EE 

riBB 

(jEB 



where N'f^ stands for the noise contribution to the 
pseudo-power spectra, a priori considered as non-zero. 
Expressions for the different noise biases under the as- 
sumption of white noise as well as those for the mode- 
mode coupling matrix depend on a choice of a specific 
formalism. They are discussed in Appendices [C] and |B 2| 
respectively, and below we only briefiy summarize main 
properties of the coupling matrices referring the reader 
for all the details there. 



1. Standard formalism 

In this formalism, assuming the same window function 
is used for estimating the E and the B pseudo-multipoles, 
the different blocks are related to each other as follows, 



Pseudo-power spectra estimators: EB 
cross-correlations 



The two types of the pseudo-multipoles defined in the 
previous Section permit us to introduce three different 
pseudo-spectrum estimators. These include a standard 
estimator in which both E and B pseudo-multipoles used 
to construct pseudo-Q's are standard; a pure estimator 
in which the E and B pseudo-multipoles are pure; and a 
hybrid estimator which combines a standard estimate of 
the E pseudo-multipoles with a pure estimate of the B 
pseudo-multipoles. 

For each of the three alternatives, the pseudo-power 
spectra, Cf^ , are defined as follows 



EE 



pEB 



m 

I TDiy) T^ivh 

m 

I V rfi;(^) r(2')* -l (04) 

2{2£+ 1) ^ [ im-^im ^ ^tm ^em\ > l^^i 



2£ 



2€ 



(22) 
(23) 



EE, EE T,^BB,BB _ urdiag 
EE,BB _ ..BB,EE ^ ..off 

-LyJ-PPI -L^J-PPI 1 



MEB.EB i\ /rdiaq 
00, Mod, 



■ m: 



off 



and 



Ti /tE E ,E B Ti /tB B ,E B r) ly/rEByEE 



(25) 



(26) 



Their explicit expressions are summarized in Eqs. (B9) 
and ( [bTqI ). 

Though M^J/ is smaller than M^^^^ ^ it cannot gen- 
erally be neglected as it describes the leakage from E- 
modes to 5-modes, which has to be corrected for, in par- 
ticular if the maps to be analyzed cover a tiny amount of 
the celestial sphere (see for example Fig. 9 of [25 ). As 
shown in App. |B 2[ the M^[,^^^ matrix is equal to zero, 
emphasizing the fact that the EB cross-spectrum is not 
coupled to the EE and BB auto-spectra. The main rea- 
son for this is that whenever the window functions con- 
tains only a £^-component, as is precisely the case in the 
standard formalism, E pseudo-multipoles are evenly cou- 
pled to E multipoles and oddly coupled to B multipoles 
while the B pseudo-multipoles are oddly coupled to E 
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multipoles and evenly coupled to B multipoles. As a con- 
sequence, the M^[,^^^ block involves product of the even 
and odd coupling and thus vanishes. This demonstrates 
that the EB cross-spectrum can therefore be treated in- 
dependently of the EE and BB spectra in this formalism. 



2. Pure formalism 

In the pure formalism, the different mode-mode cou- 
pling blocks satisfy the same relations as given in Eqs. 
(25) and (26). Nevertheless their respective expressions 



are qualitatively and quan titati vely different in both 
these cases, see Eqs. (Bll) and (B12). First of all, the 



pure estimation ensures the M^// block be small enough 
to substantially reduce amount of E-io-B leakage even in 
the cases of small-sky coverage (e.g., see Fig. 9 of [25 ). 
Second, unlike the standard case, the four blocks M^[,^^^ ^ 
which couple the EB cross-spectrum to the EE and BB 
autospectra, may not vanish. These blocks involve two 
types of couplings: i) odd-even couplings corresponding 
to the EE and BB power spectra of the window func- 
tions and ii) even-even and odd-odd couplings via the 
cross-correlation of the £^-component of the window func- 
tions with their 5-component. Though the first type of 
couplings is identically equal to zero, the second class of 
couplings does not vanish if the spin-1 and spin-2 win- 
dows have a non- vanishing 5-component. 



3. Hybrid formalism 

The case of the hybrid formalism is more intricate 
as it involves the cross-correlation of the pure pseudo- 
multipoles with the standard pseudo-multipoles. De- 
noting by X, and either E or 5, we find that 
the {EE, XX') blocks are equal to the {EE, XX') 
blocks as computed in the standard formalism while the 
{BB,XX') blocks are given by the {BB,XX') blocks 
as derived in the pure formalism. The expression for 
the {EB.XX') blocks is given by Eq. (|B13|) and in- 



volves a cross-product of the matrices Kpi?^^^^^ and 



jr^a6at//uj J ^j^^ ^^^^ ^£ ^j^^ pure formalism, the 



jdiag/off 

{EB, XX') blocks generally do not vanish as they con- 
tain even-even and odd-odd couplings corresponding to 
the cross-correlation of the £^-component of the spin-0 
window with the 5-components of the spin-1 and spin-2 
window functions. 



C. Pseudo-power spectra estimators: TB 
cross-correlations 

The formalism developed in the previous Section is eas- 
ily adaptable to the second odd-parity power spectrum, 
Cf^. Because of the E-to-B (and B-to-E) leakage, the 
TB spectrum has to be analyzed with the TE in a unified 



framework. However, cut- sky effects do not introduce 
any leakages from temperature to Q or U Stokes param- 
eter. As a consequence, their is no leakage from TT, EE, 
BB and EB spectra into the TB and TE pseudo-spectra 
and the TB spectrum can be treated independently of the 
EE,BB,EB. 

The two pseudo-spectra are defined as follows 



TE _ 



2(2^ + 1) 
1 



2(2^+1) ^ 



TprnB 



Tp^B, 



{y) 



(27) 
(28) 



from which an unbiased estimator is built by inverting 
the usual mode-mode coupling matrix, i.e. 



The explicit expressions of the mode-mode coupling ma- 
trix for the three formalisms are provided in App. |B 3| 
We observe that in the standard and pure formalisms. 



j^TE,TE 



mI^^^^ and Ml^^^^ 



-Mjpf'^^ as long 



as the same window functions are used to compute the 
E and B multipoles. However, the off-diagonal blocks, 
Mj^f'^^ are equal to zero in the standard formalism but 
may not vanish in the pure one if the spin-1 and spin- 
2 windows have a non-zero 5-component. The funda- 
mental reason is that pseudo-a^ are coupled only to 
the product of the £^-component of the window func- 
tions applied to the maps times the aj^ via an even 
coupling while in the standard formalism, £^(5)-pseudo- 
multipoles are coupled to 5(£^)-multipoles via an odd 
coupling. However, the pure formalism introduces an 
additional even coupling of the ^(5)-pseudo-multipoles 
with the 5(£^)-multipoles via the 5-component of the 
non-zero-spin window functions, which translates into 
non- vanishing off-diagonal blocks. In the hybrid formal- 
ism, the Mjpf'^^ and Mj^f'^^ blocks are formally equal 
to the ones computed in the standard formalism while 
the Mjpf'^^ and Mj^f'^^ blocks are equal to the ones 
as computed in the pure formalism. 



IV. SKY APODIZATION AND PIXELIZATION 
EFFECTS 

A. Sky apodization 

Applicability of the pure pseudo-spectrum estimator 
depends on our ability to compute sky apodization with 
the appropriate boundary conditions. Different functions 
have been proposed. They range from ones derived as a 
result of an optimization procedure to ones based on some 
analytic expressions. Their relative merit has been ex- 
tensively discussed in [25] to show that the pixel-domain 
optimization scheme proposed in [24] leads to the lowest 
variance on the power spectrum estimation. 
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The underlying strategy to derive such an optimized 
sky apodization is to search for window functions, 
which make the pure pseudo-C^ as close as possible to the 
optimal, quadratic power-spectrum estimator (see Sec- 
tion V of [24 ). An optimized weighting scheme there- 
fore involves a specific sky apodization for each £-band 
where the power spectrum is to be estimated, accord- 
ing to the signal and noise power aliasing in each band. 
Such an optimization procedure assumes that noise and 
signal are known. For the noise and £^-modes sufficiently 
precise assumptions are usually available, e.g., £^-mode 
power spectrum can be recovered precisely enough with- 
out any optimization. Though this is not the case with 
the 5-modes, it can be shown [25] that the assumptions 
about the 5-modes power spectrum very weakly affect 
the results of the power-spectrum estimation using the 
optimized sky apodization. This is due to the fact that 
the optimization process foremost attempts to reduce the 
E-io-B leakage and then the noise variance and thus no 
precise 5-mode knowledge is necessary. We note that 
in this approach as the derivative relationship is relaxed 
during the computation, the 5-component of the spin-1 
and spin-2 windows becomes non-zero. Therefore, the 
{EB^ XX) and (XX, EB) blocks of the mode-mode cou- 
pling will not vanish. Nevertheless, numerical studies 
have shown that the 5-component of their spin-1 and 
spin-2 windows is very small. As shown in App. |B 2[ only 
the cross-spectra of the £^-component and 5-component 
of the spin- weighted window functions are involved in the 
{EB^ XX) and (XX, EB) while the autospectra of both 
E- and 5-component of those window functions enter 
in the other blocks. As a consequence, we expect the 
{EB^ XX) and (XX, EB) blocks of the mode-mode cou- 
pling matrix to be much smaller than the other blocks. 

Other options for the apodizations considered to 
date [2B^ include i) analytic window functions, and ii) 
harmonic-based optimized window functions. In these 
cases, the spin-1 and spin-2 windows are computed by 
taking numerically the derivative of the spin-0 window. 
As a consequence, none of the spin-5 sky apodizations 
have a ^-component and the {EB^ XX) and (XX, EB) 
blocks of the mode-mode coupling are expected to vanish 
within numerical precision in both the pure and the hy- 
brid formalism. Hence, in these cases the EB spectrum 
is effectively decoupled from EE and BB^ as in the stan- 
dard formalism, but in a contrast with the pixel-domain 
apodization cases. 

Similar arguments apply to the case of mode-mode 
coupling for the TB cross-spectrum. In the pure 
and hybrid formalism, the potential cancellation of the 
off-diagonal blocks depends on the sky apodizations 
used in the numerical computation of the pure pseudo- 
multipoles. Those blocks are expected to be zero for 
either the analytic or the harmonic-based optimized win- 
dow functions but may become non-zero for the pixel- 
based optimized sky apodizations. From a quantitative 
point of view, the off-diagonal blocks share common fea- 
tures with the {XX, EB) and {EB.XX) blocks. They 



(TEJE) block in the standard formalism 
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(TEJE) block in the pure formalism 
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FIG. 1. The red-dashed curves corresponds to the i = 300 
column of the (TE, TE) coupling kernels computed in the 
standard (left) and pure (right) formalisms. The standard 
{TE, TE) block equals the hybrid {TE, TE) block and stan- 
dard {TB, TB) block. The overlapping black curve shows the 
TE obtained from a map containing TE and TB correlations 
dX I — 300 only. The agreement between the red and black 
curves underlines the correcteness of the numerical computa- 
tion of the different the mode-mode coupling matrices. 



in particular depend only on the cross-spectra of the E- 
component and 5-component of window functions and 
we therefore expect those off-diagonal blocks to be much 
smaller than the diagonal ones if pixel-based optimized 
sky apodizations are used. 



B. Pixelization impact on the mode- mode coupling 
matrices 

To test our implementation of the mode-mode coupling 
matrices, we compare it to the pseudo-spectrum obtained 
for an initial power spectrum with non-zero power only 
at a single angular frequency, i.e., ex For example, 

the TE pseudo-spectrum derived from an initial power 
spectrum such as Cj^ = Si^i^ and Cj^ = is equal to 

TE TE 

the elements of the £o-th column of ^, ' . Similarly, 
the TB pseudo-spectrum for such a process is given by 

TB TE 

the ^o-th column of the ^, ' block. Such a calcu- 
lation has also an additional advantage as the resulting 
pseudo-spectra include other effects, which could poten- 
tially lead to non-zero couplings, for instance, such as 
possible pixelization effects, and which are not taken into 
account in the computation of the mode-mode coupling 
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FIG. 2. The red-dashed curves corresponds to the i — 300 column of the {EB, EB) coupUng kernel for the three formalisms: 
standard, pure and hybrid formalisms (from left to right). The overlapping black curve shows the EB pseudo-spectrum 
obtained from a map containing EB correlation at £ — 300 only. The agreement between the red and black curves underlines 
the correcteness of the numerical computation of the different the mode-mode coupling matrices. 



matrices. This test has already been applied in [25 to 
the {EE, EE), {BE, BE) and {EE, BE) blocks to show 
that most of the pixelization effects are carefully taken 
into account and we will only present our results concern- 
ing the remaining blocks of the full mode-mode coupling 
kernels. 

An example of such tests is shown on Fig. [l]for the case 
of the {TE,TE), {TB,TB) and in Fig. [|for {EB,EB) 
blocks in the three proposed formalism (standard, pure 
and hybrid formalisms). The sky coverage adopted here 
is a spherical cap with a radius of ~11.3 degrees. The 
sky apodizations used are the analytic windows with a 30 
arcminutes apodization length to be applied to temper- 
ature and polarized standard pseudo-a^^'s (see [25] for 
details) and the pixel-based optimized windows for polar- 
ized pure pseudo-a^^'s. (Those peculiar sky apodizations 
are shown on Fig. |3]) In this peculiar case, a first sky 
apodization is identically applied to compute standard E- 
and 5-pseudo-multipoles and a second one is identically 
applied to compute pure £^-and 5-pseudo-multipoles. As 
already mentioned in the previous section, with such a 
weighting, the {TE, TE) block in the standard and hy- 
brid formalisms equals the {TB, TB) block in the stan- 
dard formalism, and the {TE, TE) block in the pure for- 
malism equals the {TB, TB) block in the pure and hy- 
brid formalisms. Therefore, only the {TE, TE) block in 
the standard approach and in the pure approach is de- 
picted in Fig. [l] However, the {EB, EB) block system- 
atically differs from one formalism to another and such 
a block is displayed in Fig. [2] for the three techniques. 
In such figures, the red-dashed curves correspond to the 
£ = 300 column of the computed mode-mode coupling 
while the overlapping black curves stand for TE{TB) 
and EB pseudo-spectra as derived from a map contain- 
ing correlated T-, B- and £^-modes at ^ = 300 only. The 
overall agreement of red and black curves pins down the 
right computation of the mode-mode coupling matrices. 
The agreement is nevertheless not perfect for two reasons. 
First, the equality between the column of the mode- mode 
coupling and the pseudo-spectrum from a single mode 



process is exact on average. However, the black curves 
displayed on Figs, [l] and [2] are obtained from one real- 
ization only, explaining the scattering of the black curves 
around the red ones. Second, we can not exclude the 
influence of residual pixel effects. 

As already explained, the {TB,TE) and {EB,PP) 
blocks are non-vanishing (though expected to be small) 
in the peculiar case of the pure and hybrid formalisms 
if and only if the spin-1 and spin-2 sky apodizations 
have a non-zero 5-component. This is the case of 
pixel-based optimized sky apodizations. However, the 
above described tests shows that within the HEALPix 
scheme, and irrespectively of any formalisms, unmod- 
elled pixel effects induce some non- vanishing {TB, TE) 
and {EB,PP) blocks with an amplitude higher than 
the cut-sky effects expected in the the pure and hybrid 
formalismq^ However, the pixel-induced non- vanishing 
{TB, TE) and {EB, PP) blocks in the HEALPix scheme 
remains much smaller than the cut-sky (and modelled) 
{TE,TE) and {EB,EB) blocks. We will therefore ne- 
glect such effects and systematically set the {TB, TE) 
and {EB, PP) blocks to zero in the following. 

As the pixel-induced leakages are not corrected, the fi- 
nal estimate of all the angular power spectra but the TT 
one is a priori biased as long as the TB and EB correla- 
tions do not vanish. Fortunately, such cross-correlations 
are expected to be zero in the standard cosmological 
paradigm. Moreover, we have tested that with non- 
vanishing TB and EB input power spectra at a level of 
l{l^l)Cj^ /{2ii) - 0.2 iiK^ and £(^ + l)Cf ^/(27r) - 
0.004 iiK'^ at £ = 2- 1000 (as motivated by the potential 
impact of an helical, primordial magnetic field [HI [15]), 



^ As is already shown in [25] (Sec. IV. A. 5) for the specific case of 
EE and BB spectra , such unmodelled pixelization effects are 
significantly reduced by using ECP-like pixelization schemes (for 
which the number of azimuthal grid points remains unchanged 
for each constant zenithal ring and therefore providing a nearly 
exact quadrature over the polar angle). 
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the EE^ BB, TE, TB and EB angular power spectra es- 
timated on 1% of the sky with inhomogeneous noise (see 
Fig. |8] for the noise distribution) are indeed not biased 
though residual EB ^ EE, EB ^ BB and TB ^ TE, 
leakages are not corrected for [4Qj . 



V. APPLICATIONS TO SMALL-SKY 
EXPERIMENTS: POWER SPECTRA 
UNCERTAINTY 

In this Section we discuss the performance and appli- 
cability of the above defined odd-parity power spectra 
estimators in the context of vanishing true TB and EB 
cross-correlations, i.e., considered as a null test, poten- 
tially used for search of astrophysical and/or instrumen- 
tal systematic effects. The uncertainty of each estimator 
is quantified by computing the power spectrum variance 
thanks to Monte-Carlo simulations. Our fiducial model 
for the input T and E CMB signal is given by the cosmo- 
logical parameters as constrained by the WMAP 7-year 
data. The input B-mode is composed of a primordial 
component with a tensor-to-scalar ratio r = 0.05 and 
a secondary component induced by lensin^ The two 
input odd-parity power spectra are set to zero (as is pre- 
dicted in a parity invariant cosmological scenario). 

Our starting point is an idealized mock survey con- 
sidered below covering a circular patch of ~1% of the 
sky area with a level of homogeneous noise set to 5.75 
/iK-arcmin for the three Stokes parameters. In order to 
get closer to realistic sky observations, two possible ex- 
tensions of this ideal mock survey are implemented. We 
first consider the effect of more intricate boundaries us- 



ing a square patch with holes due to e.g. point-sources 
removal, and later an implementation of an inhomoge- 
neous sky survey as expected for balloon-borne experi- 
ment. The sky coverage is always taken to be ~1% and 
the RMS of the noise in the inhomogeneous case is set 
equal to 5.75 /iK-arcmin. The values of both these pa- 
rameters are motivated by the ebex experiment [32 . In 
addition, in all cases we assume that noise is uncorre- 
lated from pixel to pixel and from one Stokes parameter 
to another, i.e., 

{Ns{i)Ns^{j))=4m^j6s,s'. (29) 

As a consequence, the TE, TB and EB noise biases van- 
ish (see Appendix [C]) . 

Hereafter, the variance is estimated as the standard 
deviation of 500 MC simulations. The input {CMB + 
noise} maps are built using the HEALPix pixelization 
scheme at a resolution of N^ide = 512 corresponding to 
a pixel size of ~ 7 arcmin. The estimated power spectra 
are binned with a band-width of A£ = 40 and with the 
lowest bin starting at ^ = 20. The pseudo-a^^ for tem- 
perature are estimated applying the arc cosine apodiza- 
tion [25] with an apodization length of 30 arcminutes. 
The same apodization is used while calculating the stan- 
dard pseudo-multipoles for polarization from the Q and 
U maps. For the pure-pseudo-multipoles of E- and B- 
modes computation we use a pixel-based, optimized sky 
apodization to the Q, U maps. 

As a reference against which to compare the variances 
obtained from the MC simulations we use a theoretical - 
and optimistic- Fisher estimate of the variance, which for 
a cross-spectrum of two sets of (/, Q, U) maps, labelled 
by (A) and (5), read^ 
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Here, A 



^x^A)y{B) 



denotes the variance 



of a cross-spectrum of -^(a) and Y^b) where X and Y 



Our convention for r follows the WMAP convention: r = 
^T(ko)/^s{ko) with 7^s(T)5 the primordial scalar(tensor) power 
spectrum and ko = 0.002 Mpc"-*^ the pivot scale. 
We point out that in [25], a factor 1/2 is missing in front of the 
last term (pure noise term) of their /g^y -formulae. Their equation 
(32) for cross-spectrum should read 



ACf 



1 

sky 



\ (2^+1) 



B2 



1^4 

2Bj 



stand for either T, E oi B as derived from a map either 
A oi B respectively, A/'pix = 12N'^:^^^ is the total number 
of pixels, and 5£,(a) stands for a resolution of the map 
A. Given our assumptions about the noise, Eq. (29), 

2 2 r 

(7 V \^ — O' V V ^ Ox Y 



where g\ 



is related to the noise per pixel crp.^(z) via a 
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,sEi^pix(0- The effective sky fraction, is 
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defined as 



Apodized mask 



Spin — Window 
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„X(A)y(B) 

^sky 



i=l 



N ■ 

-L^pix 2 



(31) 



i=l 



where ^^^ ^""^ = W^^^^ x W'^''\ is a cross-product of 
tiie sky apodizations applied to tiie maps A and B from 
wiiicii tiie X- and F-modes are estimated. As ttie sky 
fraction depends on tiie ciiosen sky apodizations, wiiicii 
in general will depend on the adopted formalism, on the 
map itself, and the considered ^-bin, the resulting Fisher 
estimate of the error bars defined in Eq. (30), will also 
defined on all these factors. As a unique -formalism 
independent- benchmark, we will therefore use the Fisher 
variance computed for the maximal value of /sky for each 
considered map corresponding to setting the sky apodiza- 
tion equal to the binary mask. This theoretical uncer- 
tainties therefore corresponds to the lowest Fisher esti- 
mate one may expect given a peculiar sky fraction, not 
taking into account specific apodization used in the cal- 
culations of the corresponding pseudo-spectra. 

The entire set of operations needed to compute the 
estimates, i.e., (1) the mode-mode coupling matrices, 
(2) the six pseudo-cross-spectra, and (3) the six angu- 
lar cross-power spectra estimates from the noise-debiased 
pseudo-C^, using CMB maps and sky apodizations as in- 
puts, are numerically implemented in the X^pure code. 
The code is an extension of the Xspect and Xpol codes 
[28] and is based on the (pure)s^hat library -an effi- 
cient massively parallel implementation of spin- weighted 
spherical harmonic transforms [36H39] . This code can 
also be used in an MC setting allowing the user to sim- 
ulate, and subsequently analyze, random realizations of 
CMB intensity and polarization maps from the six angu- 
lar power spectra as an input model. 



A. Results for an ideal mock survey 

The idealized mock survey consists of a circular patch 
with a radius of ~11.3 deg., centered at the equator to 
reduce pixel effects (see section IV. A. 5 of [25 ), and with 
homogeneous noise at a level of 5.75 /iK-arcmin. Such 
a survey is ideal due to both the simplifying assumption 
of noise homogeneity and its simple shape. Indeed, a 
circular patch exhibits the shortest perimeter for a fixed 
value of the sky coverage. As ^-to-^' and E-io-B mix- 
ing appears as the result of boundary effects, such shape 
permits to minimize different mode-mode couplings for 
a given sky coverage. In addition, we do not allow here 
for any internal boundaries, e.g., due to masked point 
sources. The analytic sky apodization -to be applied 
to intensity map and polarization maps for standard 
pseudo-a^^'s-, and the three spin- weighted sky apodiza- 
tion, optimized for the second bin, I G (60, 100) -and 
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FIG. 3. Analytic sky apodization (to be applied to inten- 
sity map and polarization map for standard pseudo-a^m's es- 
timate) and spin-0, spin-1 and spin-2 sky apodizations, op- 
timized for the second bin, I G (60, 100) (to be applied to 
the polarization maps for pure pseudo-a^rn's estimate), in the 
case of the idealized mock survey. (Only the real parts of the 
spin-1 and spin-2 windows are shown.) 



to be applied to the polarization maps for pure pseudo- 
<^^m's-, are displayed on Fig. |3] 

For this ideal mock survey, the fraction of the sky 
is /sky — 0.96%. For the TT power spectrum (all for- 
malisms), the EE and TE power spectra (standard and 
hybrid formalisms), and the BB^ TB and EB (standard 
formalism), the effective sky coverage is constant and 
equal to 0.92%. For the BB power spectrum (pure and 
hybrid formalism), and the EE and EB power spectra 
(pure formalism), it varies from 0.64% at large angular 
scales (20 < ^ < 60) up to 0.95% at small angular scales 
(980 < I < 1020). Finally, for the TB power spectrum 
(pure and hybrid formalisms), the TE power spectrum 
(pure formalism), and the EB power spectrum (hybrid 
formalism), the sky coverage ranges from 0.78% at large 
angular scales (20 < ^ < 60) up to 0.93% at small angular 
scales (980 <i< 1020). 

The MC averaged estimate of the six power spectra 
{i.e. TT, EE, BB, TE, TB and EB spectra), esti- 
mated using the three types of formalism, are displayed 
on Fig. [4] alongside the theoretical input C^'s. The per- 
fect agreement between the estimated and the theoretical 
input shows that each of the formalism indeed leads to an 
unbiased estimate of the angular power spectra. As such 
estimators are built to be unbiased, this essentially con- 
firm the correct numerical implementation of the three 
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FIG. 4. MC averaged estimation of the six angular power spectra (colored curves) alongside the theoretical input C^'s (black 
curves) for the case of an ideal mock survey consisting of a sphercial cap with homogeneous noise (see Fig. [s]). The perfect 
agreement between the estimated and the theoretical C^'s proves that each of the three formalisms leads to an unbiased 
estimation. 



formalisms in the X^PURE code. 

The statistical uncertainties of the estimated angular 
power spectra obtained via MC simulations are displayed 
in Fig. [5] alongside the Fisher estimates of the variance 
computed as explained in the previous Section, Eq. 
(30). We emphasize that except for TT, all the studied 



formalisms are theoretically different from each other, as 
for instance they lead to different mode-mode coupling 
matrices, and thus could result in different levels of the 
estimates uncertainties. Nevertheless, as the (TE^TB) 
block is typically negligible, we can expect that the TE 
estimate will be identical in the standard and hybrid 
formalisms, and that the TB estimates in the pure and 
hybrid formalism will coincide. For the EE^ BB and 
EB spectra, the mode-mode coupling matrix changes 
from one formalism to another, even if {EE, EB) and 
{BB, EB) blocks are neglected. We consequently expect 
differences -though the hybrid and standard estimates 
should be close in the case of the EE spectrum and 
the hybrid and the pure estimates of the BB spectrum 
should also be rather similar as well. Some other 
expectations are as follows: 

• TT-spectrum: For such an angular power spectrum, 
the three formalism are equivalent and we expect to 
reach the same level of uncertainties for all the formal- 
ism. This is clearly confirmed by the results exhibited in 
the top-left panel of Fig. [Sj 

• EE-spectrum: As already suggested in [23 and illus- 
trated in [26], the standard estimator leads to slightly 
smaller error bars than the pure estimator. The under- 
lying reason is that, though B-mode leaks into T^-mode 



in the standard approach (thus increasing the sampling 
variance), this effect is small and the information lost by 
removing ambiguous mode in the pure approach leads to 
a higher increase of the sampling variance. In addition, 
the hybrid estimator, being based on the standard 
pseudo-multipoles for T^-modes, performs roughly the 
same than the standard estimators. The gain obtained 
by using the standard or hybrid estimators is clearly 
seen on the top-middle panel of Fig. [5] at large angular 
scales. 

• BB-spectrum: The case of BB power spectrum is 
precisely the opposite of the EE one. The main source 
of uncertainties is the extra-sampling variance due to 
£^-modes leaking into 5-modes. As clearly illustrated on 
the top-right panel of Fig. [5j both the pure and hybrid 
estimators perform the same, leading to uncertainties 
on par of the most optimistic Fisher estimate, while the 
standard estimator leads to uncertainties greater than 
Cf^ for I smaller than - 300. 

• TE-spectrum: For the very same reason that standard 
and hybrid estimator leads to smaller uncertainties 
than pure estimator for EE spectrum, one expect 
those two approaches to also perform the best for TE 
cross-spectrum. Results displayed on Fig. 5] (bottom- left 
panel) show significant improvement by using the 
standard or hybrid estimator as compared to the pure 
ones especially at large angular scales. 

• TB-spectrum: This case is conceptually identical 
to BB power spectrum. The lowest uncertainties are 
achieved by removing T^- modes leaking into B: pure or 
hybrid formalism are then preferred to get an estimate 
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FIG. 5. Error bars on the reconstructed angular power spectra for each of the three formahsms (colored curves) alongside the 
naive (binned) Fisher estimate of such uncertainties. The assumed sky survey is the ideal case consisting of a circular patch 
with homogeneous noise (see Fig. [sl. The error bars are obtained as the standard deviation over 500 MC simulations. 



of Cj^ as accurate as possible with respect to statistical 
error bars, as exhibited on the bottom- middle panel of 
Fig.|5l 

• EB- spectrum: Standard estimators, pure estimators 
and hybrid estimators perform the best respectively 
for £^-pseudo-multipoles only, for 5-pseudo-multipoles 
only and for both E- and 5-pseudo-multipoles re- 
spectively. Our results show that the gain in using 
pure B-pseudo-mulipoles is more significant than the 
gain in using standard £^-pseudo-multipoles. One 
then expect the following hierarchy for the error bars: 
Standard > ^f^ure > iybrid' Though Uncertain- 
ties from the standard estimators are clearly higher (see 
bottom-right panel of Fig. [5|, it appears that the pure 
and hybrid estimators perform the same. 

From these results, we could conclude that the hybrid 
estimator systematically leads to the lowest uncertainties 
for the six angular power spectra; an expected conclusion 
as the hybrid formalism is based on the best pseudo- 
multipoles computation for both E- and 5-modes. 



B. Towards a realistic survey: effect of holes and 
sharp boundaries 
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FIG. 6. Same as Fig. [3] for the case of a mock survey with 
sharpenned boundaries and holes. 



In this Section we consider a more intricate observed 
patch shape consisting of a square patch with internal 
holes mimicking effects of point source masking [24 . The 
specific patch used hereafter is depicted in Fig. [6j 

The fraction of the sky covered by this patch is /sky — 



0.99%. For the TT power spectrum (all formalisms), the 
EE and TE power spectra (standard and hybrid for- 
malisms), and the BB^ TB and EB (standard formal- 
ism), the effective sky coverage is constant and equal 
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FIG. 7. Error bars on the reconstructed angular power spectra for each of the three formaUsms (colored curves) alongside 
the naive (binned) Fisher estimate of such uncertainties. The assumed sky survey is a square patch filled with holes and with 
homogeneous noise (see Fig. [6|.The error bars are obtained as the standard deviation over 500 MC simulations. 



to 0.92%. For the BB power spectrum (pure and hy- 
brid formalisDi) , and the EE and EB power spectra 
(pure formalism), it varies from 0.43% at large angular 
scales (20 < i < 60) up to 0.96% at small angular scales 
(980 < £ < 1020). Finally, for the TB power spectrum 
(pure and hybrid formalisms), the TE power spectrum 
(pure formalism), and the EB power spectrum (hybrid 
formalism), the sky coverage ranges from 0.78% at large 
angular scales (20 < i < 60) up to 0.93% at small angular 
scales (980 <£< 1020). 

We found MC averaged power spectra in perfect agree- 
ment with the theoretical input Ci showing that power 
spectra estimators for the three formalisms are unbiased 
as expected. 

The error bars on the estimates of the six angular 
power spectra are shown in Fig. (T] As already shown 
for the case of the circular patch, the standard or hy- 
brid estimators gives smaller variance for the EE and 
TE power spectra, while the pure or hybrid formalisms 
are to be adopted for the BB^ TB and EB spectra. One 
should however expect some differences: since sharpened 
and internal boundaries enhance the E-io-B and B-io- 
E leakages {i.e. the total power contained in ambiguous 
modes is increased), the relative merit of each formalism 
is supposed to be quantitatively strengthened. Though 
not obvious for BB and TB spectra (for which the gain 
in using the pure or the hybrid estimators is already dra- 
matic for simple boundaries), it is now clear at large an- 
gular scales {£ < 100) that the error bars on EE and 
TE angular power spectra (see top-middle and bottom- 
left panels of Fig. [t]) are lower by using the standard or 



the hybrid formalism than by using the pure approach. 
Moreover, a closer inspection of the EB power spectra 
seen in the bottom-right of Fig. [7| suggests that, at large 
angular scales, £ < 150, the hybrid approach performs 
slightly better than the pure approach, as expected from 
theoretical considerations. 



C. Towards a realistic survey: effect of 
inhomogeneous sky sampling 

For a realistic small-scales experiment, the observed 
patch is expected to be rather irregular with the density 
of observations per sky area strongly varying from one 
pixel to another. This will unavoidably lead to signifi- 
cant noise inhomogeneity of the reconstructed sky maps. 
We therefore explore the performance of the TB and EB 
pseudo-C^'s estimators in the case of an observation mim- 
icking a long-duration CMB balloon-borne experiment 
such as EBEX [32,. The resulting sky sampling, i.e. a 
number of observations per pixel (Fig. [8|, has already 
been used in [2F to assess the performance of the pure 
B-mode estimators. It ranges from 10^ at the edges of 
the patch to more than 2 x 10^ in its core after we cut 
the noisiest pixels in order to ensure convergency of the 
optimized sky apodization procedure (see [25]). Finally, 
we set the noise per sample such as the average noise 
level for the full patch equals the noise level used in the 
previously-studied homogeneous cases. 

The retained sky fraction is equal to /sky — 0.89%. 
For the TT power spectrum (all formalisms), the EE 
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Complete distribution of sky observations 
1 



Cut distribution of sky observations 



Cut distribution of si^y observations 




FIG. 8. A distribution of sky observation for an observational strategy mimicking the ebex experiment. The complete 
distribution is shown in linear scale (left panel) for which the density of samples ranges from 1 to ^ 2 x 10^. The noisiest 
pixels, which lie in the outer part of the observed sky, have been trimmed in order to compute the optimized sky apodizations. 
This "well-observed" part is displayed in linear scale (middle panel) and in logarithmic scale (right panel). The density of 
observations for the conserved pixels ranges from 10^ to 2 x 10^. 



and TE power spectra (standard and hybrid formalisms), 
and the BB^ TB and EB (standard formalism, the ef- 
fective sky coverage is constant and equal to 0.85%. For 
the BB power spectrum (pure and hybrid formalism), 
and the EE and EB power spectra (pure formalism), it 
varies from 0.42% at large angular scales (20 < f < 60) 
down to 0.29% at small angular scales (980 <l< 1020). 
Finally, for the TB power spectrum (pure and hybrid 
formalisms), the TE power spectrum (pure formalism), 
and the EB power spectrum (hybrid formalism), the 
sky coverage ranges from 0.6% at large angular scales 
(20 < I < 60) down to 0.46% at small angular scales 
(980 < i < 1020). We note that for this inhomoge- 
neous noise case, the effective sky coverage decreases 
for smaller scales as the optimized sky apodization con- 
verges towards inverse- noise weighting, which is "natu- 
rally" apodized as a consequence of the scanning strat- 
egy. This difference of the apodization lengths can be 
easily seen by comparing the left panel of Fig. [Sj show- 
ing the inverse noise weighting, with the right top panel 
of Fig. [9) depicting one of the optimized, large scale 
apodizations. 

Once again, we first checked that each formalism leads 
to unbiased estimates for the six power spectra (TT, EE^ 
BB, TE, TB, EB). 

The uncertainties of the estimated angular power spec- 
tra are shown in Fig. 10 alongside the Fisher estimates 
Overall the hybrid formalism leads to the smallest error 
bars and allowing for the inhomogeneous noise does not 
change at a qualitative level, the main conclusion drawn 
earlier based on the homogeneous noise cases. A similar 



Apodized mask 



Spin — Window 
e=(60,ioo) 



^ We note here in passing that the difference of the BB pure spec- 
trum shown here with that shown in Fig. 24 of (25| is due to a 
different noise level assumed erroneously in this earlier work. 




FIG. 9. Same as Fig. [3] for the case of a realistic mock survey 
with inhomogeneous noise as displayed in Fig. [S] 



conclusion applies to the relative merits of each formal- 
ism. On the quantitative level, the noise inhomogeneities 
seem to enhance the differences between the techniques, 
as, for example, can be seen clearly in the case of the EE 
and TE power spectra in Fig. 



10 The larger uncertain- 



ties seen in these cases for the pure estimator are due to 
the fact that sky apodizations employed in the pure esti- 
mation of ^-modes are a priori not optimized to account 
for the noise variation. 
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20 100 1000 20 100 1000 20 100 1000 

Multipole t Multipole I Multipole t 



FIG. 10. Error bars on the reconstructed angular power spectra for each of the three formahsms (colored curves) alongside the 
naive (binned) Fisher estimate of such uncertainties. The assumed sky survey is shown in Fig. [9] with inhomogeneous noise as 
displayed on Fig. ^.The error bars are obtained as the standard deviation over 500 MC simulations. 



VI. CONCLUSIONS 



We have presented three unbiased alternative ap- 
proaches to estimate the six angular power spectra associ- 
ated to CMB temperature and polarization anisotropics. 
Those techniques are based on the so-called pseudo- 
C^'s estimators and generalize to the full set of angular 
power spectra the approaches already proposed in [23]- 
[25l l27l l28j . Each technique differs from the other in the 
way E- and 5-pseudo-multipoles are computed from the 
maps of Stokes parameters. The three approaches are 
equivalent from the perspective of the TT power spec- 
trum. The first approach, dubbed standard^ estimates 
polarization multipoles by projecting the masked Q and 
U maps on the standard spherical harmonics of E- and 
B-type. This technique does not lose any information 
contained in the observed maps but its variance suffers 
from the E-io-B and B-io-E leakages. The second ap- 
proach, denoted pure^ makes use of the pure-£^, pure- 
B and ambiguous modes decomposition [35l to project 
the Q and U maps on a E- and 5-basis free from any 
leakages. This second technique therefore corrects for E- 
to-B and B-io-E leakages on any single realization but 
loses some part of the information by removing ambigu- 
ous modes. Finally, in the third approach, called hybrid^ 
^-pseudo-multipoles are numerically computed using the 
standard approach while 5-pseudo-multipoles are com- 
puted using the pure technique. For the specific case of 
CMB polarized anisotropics, ^-modes are much greater 
than 5-modes. From theoretical arguments and previous 
numerical experiments, statistical uncertainties on the 



£^-pseudo-multipoles are therefore expected to be mainly 
increased by the loss of information while statistical un- 
certainties on the 5-pseudo-multipoles are expected to 
be mainly increased by E-to-B leakage. As a result, the 
hybrid technique should provide the smallest error bars 
at least as long as the B-mode power is not larger than 
that of the £^-modes. 

We have developed and implemented those three ap- 
proaches by providing the set of expressions to efficently 
compute the different mode-mode coupling kernels for 
each of the formalisms, as well as the noise bias for the 
specific case of uncorrelated noise from pixel to pixel. 
The overall consistency of those computations has been 
demonstrated via numerical experiments. 

Finally, we have assessed the relative performances of 
the three proposed pseudo-Q's estimators for the case 
of small-scales experiments, assuming three different ex- 
perimental configurations. In all these cases and for all 
6 angular power spectra, the hybrid estimator has been 
found to result in the smallest statistical uncertainties, 
confirming the theoretical expectations, while the other 
estimators can produce similar uncertainties in some spe- 
cific cases. For instance, i) the standard and hybrid for- 
malisms performs equally well for estimating the EE and 
TE power spectra, ii) the pure and hybrid formalisms 
performs equally well to estimate the BB and TB and 
iii) the hybrid technique performs the best (at least at 
large angular scales) for the EB power spectra, though 
the pure approach is on par with the error bars achieved 
with the hybrid estimator. 
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Appendix A: Gaunt integrals and Wigner-3j 

We provide in this appendix some general definitions and properties of the Gaunt integrals and Wigner-3j symbols 
which are useful for the computation of the mode-mode coupling matrices of the pseudo-C^ estimators. 
The general definition of the Gaunt integrals is 



/ s^^m s^ytm^ s"yi"ni"dn, (Al) 

and it can be expressed by use of the Wigner-3j symbols 



niras 

^ I' m' s' m" s" 

f4n 



rims _ . / (2£ + me' + i)(2r + 1) fe e r \ f i e r , 

yi'm's'-x"m"s" - y 4^ \s s' s" )\m m! m!' )■ 

Ths is easily generalized to integrals involving complex conjugates of the spin- weighted spherical harmonics by using 
s^lm ~ {~^y'^"^ -sYi(-m)- The Wigner-3j symbols verify the following symmetry and orthogonality relations 



and 



mm' 

Two linear combinations of Wigner-3j symbols and Gaunt integrals are involved in the expression of the pseudo- 
multipoles mode-mode coupling matrices. To lighten our expressions, we defined the following quantities: 



and 



Gitl{L,L',L")= I 



{sylm^y^'m' -2+8^1" m") i (-s^^^-2^rm' ) 



= (-1) 



47r _ 

s+m 



dn, (A6) 



/ (2^ + i)(2f + i)(2r + i) / 1 e r \ (±) , „ 



with L labelling the couple (^, m) in the Gaunt integrals. From the parity relation of the Wigner-3j, it is easily shown 

that 3^2^^ = for odd values of + f ' + i") and J^Z] = for even values of + + Finally, we also define the 
following Wigner-3j symbols 



which is equal to zero for odd values of {£-{-£'■ 
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Appendix B: Mode-mode coupling matrices 

1. Mode- mode coupling for (pure) pseudo-a^rn 

In this appendix, we present the exphcit expressions of the different mode-mode couphng matrices related the 
pseudo-multipoles to the CMB multipoles. We remind that standard pseudo-a^^'s are defined by 



B 



(std) 



' ^ w 



w 



and pure pseudo- 



^'s are defined by 

-^{pure) 1 



(Q + iU) [d'^Yi^y + {Q- iU) {B'^Yi^] 
{Q + iU) {d^Ye^y - (Q - iU) {B^Yemy 



dn, 



dn. 



El 



{pure) 



i 



27V, 



',2 J47r 



(Q + iU) {d^WYemY + (Q - iU) {B^WY^^y 
(Q + iU) {d^WYemY -{Q- iU) [d'WYtrr 



dn. 



(Bl) 
(B2) 

(B3) 
(B4) 



We will first provide the expressions for pure pseudo-multipoles as the standard approach can be easily deduced from 
the pure one. 

a. Pure pseudo-multipoles 

By expanding the spin-raising (spin- lowering) operators acting on the product W x Yirn and using the definition of 
the spin-weighted spherical harmonics, one obtains 



{pure) 
im 



B^ 



{pure) 



-It. 

/ E 

2 Jatt „_n 



a,Np 



(Q + iU) {d^-^wy,Yi^ + {-ly {Q - iU) {d^-^wy 
(Q + iU) {d^-^wy - (-1)^ (Q - iU) {B^-^wy .,Y, 



dn. 



dn. 



(B5) 



(B6) 



with = 1 for 5 = 0, 2 and = 2 for 5 = 1. The 5 = 2 term in the above summation over the spin index simply 
corresponds to the standard multipoles. Therefore, one can first compute the mode-mode coupling matrices for the 
pure case and then keep only the s = 2 term in the summation to recover the mode-mode coupling matrices for the 
standard estimation. This computation is done by plugging the multipolar decomposition of the (Q ± iU) fields and 
of the window function W , i.e. 



W = -Y.wl%Ye^, 

£m 

Q ±iU = - ^ {aE/m ± idB/m) ±2Y£m- 



If the derivative relationship between the different spin-weighted window function is completely statisfied, i.e. Ws = 
d^W, the multipolar decomposition of the spin-0 window function W is sufficient as 

^SL = ^i,sWim and = 0- 

However, as shown in [2F, pixelization on the sphere partially breaks this derivative relationship and the multipolar 
decomposition of each of the spin-weighted window function has to be considered as independant from the spin-0 
ones, i.e. 



im 



18 



By plugging the above multipolar decomposition into Eqs. (B5) and (B6) and using the Gaunt integrals, it is 
straightforward, though rather long, to show that 



with 



[pure) 



T^ipure) _ \ rrdiag nTT^ff 

i'm' 



Tj-diag 



-1 



off 

£m,i'm' 



-1 



We underline that the spin-dependant pseudo-multipoles, Es^tm and Bs^tm-, as defined in Eqs. (18) and (19) are 
also linked to the CMB multipoles via a mode-mode coupling 



TTi \ ^ jys,diag , • rySjOff 



with 



s,diag 



K. 



^ {-lya, [wfX,.^..G'C}s{L, L\ L") - iw[%„^„G^-},{L, L', L") 

l"m" 

^ {-lya, U%,,^,,g[-1{L, L\ L") - iwi''_le"m"Gt\iL, L', L") 



The full mode-mode coupling matrices are then related to the spin-weighted mode-mode coupling matrices by 

1 ^ 1 ^ 

-— Nn K^^dzag . TjOff _ ^ yvr T^s,off 



diag 



b. Standard pseudo-multipoles 



From the above result, we can easily deduce the mode- mode coupling matrices for standard multipoles by simply 
keeping the 5 = 2 terms in the summation over the spin index to get 



E 



(std) 



B 



(std) 
im 



I'm' 

E\ rrdiag • r^off 



i'm' 



with 



diag 



K. 



em, i'm' 



2, diag 



K 



K'^^off 

-^^im/'m' 



\Y.<rm"Gt\{L,L\L"\ 

i"ra" 

\Y.<i"m"G[-\L,L',L"). 



(B7) 



l"ni" 



We stress out that for standard estimation of the pseudo-multipoles, only the real-valued, spin-0 window function is 
used which does not contain any 5-component. 
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2. Mode-mode coupling for EB pseudo-C^ 



We remind that the pseudo-multipoles are related to the CMB multipoles via a mode-mode coupling matrix 



E 



(x) 



E 



\ i^diag 



B 



{y) 



H 



diag 
tm-1'- 



£'m' 



where (x) and {y) corresponds to either "standard" or "pure" estimation. We stress out that the above expressions 
do not formally depend on the formalism. All the formalism dependancy is encoded in the matrices /C and which 
can be set equal to K or depending on the values adopted for {x) and {y). In the following, we assumes that the 
polarization field is a gaussian and isotropic random fields 



{ciE,£niCiB,£'m') = B ,£niCiE /' m') = ^ £ ^£ 



r>EE 

^£,{xx) 

pBB 



2i- 
] 

2i- 



£ra^£ra ' 



The pseudo-power spectra are defined as follows: 

m 

J TD{y) T^iy> 

1 ^ / .^£m^£m 

2(2^ + 1) ^ 



fiEB 

^£,{xy) 



'm-^£m 



■ ^£m^£m 



Starting from the definition of the pseudo-multipoles, it is easily shown that 
1 



nEE 

^£Xxx) 



2i- 



1 



BB \ ^ 



diag 



EE 



off 

£m-.£'m' 



cf; 



EE 



off 



BB 



diag 
£m-.£'m' 



' r^BB 



2Im 



2Im 



diag 



-off-k 



:m:l"m' '^£m:l"m' 



EB 



r>EB 

^£,{xy) 



1 fp^ njdiag-k _ ^off njoff-k 1 ^E 

Y Y^£m;£'m'^£m;£'m' £m;£' m' ^ £ra]£' ra' \ ^£' 



^diag njoff-k 



Cfi 



r^EB 



EB 



2e 

+im \n 



(B8) 



£' mm' 
diag i^off-k 



£m]£' m' £m]£' m' 



^ 0/ 



BB 



Im 



jy^diag njoff-k 
'^£m-l'm' ^£m-l'm' 



EE 



}■ 



The form of the above expressions is also formalism-independant as all the specificities of a given formalism is encoded 
in the pect^/mr expressions of the /C and 1-L matrices. Without any a priori on the explicit expressions of those matrices, 
each of the three pseudo-C^'s receives contribution from the three angular power spectra. However, due to the parity 
properties of Wigner-3j symbols, it will be shown in the following of this appendix that the {EB^ PP) and (PP, EB) 
blocks are zero as long as the window functions do not contain any P-type component. 

The mode-mode coupling matrices for the pseudo-power spectra are finally computed by taking the expression of 
the /C and 1-L matrices valid for the adopted formalism and then performing the summation over the two azimuthal 
quantum numbers m and m' . This explicit computation is very similar to the temperature case [27l |33l |34]. The 
standard and pure formalism being a restricted case of the pure formalisnFj the computation of the mixing kernels 
is only detailed for this last case. The forthcoming formulae involve the J^^{l^l' ^l") quantities and we will remove 
their ^' and I" dependance to lighten the writings. 



^ We remind that the pure pseudo-multipoles are a linear combi- 
nation of three spin-dependant pseudo-multipoles (s = 0, 1, 2) 



and the standard pseudo-multipoles are given by the spin-2 term 
of the pure ones. 
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a. Standard formalism 



In the standard formalism, the mode-mode couphng matrices for the pseudo-spectra reads: 

167r 



^EE^EE _ j^BB^BB _ + 1 V- \ (E) .(+) 



i\jEE,BB 



i\/rEE,EB 
^^J- ppi 



,rBB,EE _ + 1 



(-) 



2f 



t"m' 

-'-y-'- ppi 



L(s) j(+)r_ ^(s) j(-) 



2e 



^0 ^0 ' 



EB,BB _ + 1 ^ I (£;) 



^0 ^0 



(B9) 



However, because J^+^ vanish for odd values of {£ -\- £' -\- 1") while ^ vanish for even values of + l")^ the 
off-diagonal blocks relating autospectra to the EB cross-spectrum become null 



BB,EB 



EB.EE 



EB,BB 



0. 



(BIO) 



b. Pure formalism 



For the {EE^ EE) block computation, the expression of Gaunt integrals as functions of Wigner-3j symbols is first 
plugged in the expression of Hf!^%^f to get 



..XE) t(+) (p pi pm\ _ /(-) (p pf piii\ 



We can now first use the orthonormalization properties of the Wigner-3j listed in App. [A] and secondly, the fact that 
all the terms proportional to (^^2-s^2-s') equal to zero to get 



.yrEE,EE _ + 1 v-^ 



E"' 

s=0 



w. 



( + ) 



2-sl"m"'^2-s 



E 

s=0 



a. 



Np. 



(-) 



2-s,£"m"'^2-s 



(Bll) 



We emphasize that the spin-summation runs from to 2 for the £^-component of the window functions but from to 
1 for its ^-component as only the spin-1 and spin-2 windows may have a non- vanishing 5-part. Applying the same 
reasoning for the {EE^ BB) and {EE^ EB) blocks leads to 

2 1 21 



lyjEE.BB _ + 1 



E 

2 



a, 



(-) 



2-s,£"m"'^2-s 



s=0 



J 



( + ) 



2-sl"m"'^2-s 



i^/rEE,EB 2£^-\-l 

Mpp, = — 2^ 2^ asas 



Stt 



=0 



^2 



W 



2-s,l"ni"^2-s' ,l"m" 



^2-s^2-s' [^2-s,^''m'''^2-s',^''m"J J * 

The other blocks are easily deduced from the above three ones as 



(B12) 



MB B ,B B -j\yrEE,EE 
ppi — ^^J-pp' •) 



jyjBB,EB 



][ /fE E ^E B -n ^E B ,E E 



-\ /rE E ^E B -n ^E B ,B B -n ^E E ,E B 
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The expression of the Mf^f'^^ block shows that in the pure formahsm, the EB-EE couphng and the EB-BB couphng 
are not zero (unhke the standard formahsm). However, such couphngs wih vanish if the spin-1 and spin-2 window 
functions have a vanishing component of type B. 

c. Hybrid formalism 

In the hybrid formahsm, only the MfJ^'^^ blocks needs to be computed as the Mf^f'^^ blocks are equal to the 

EE P P' BB P P' BB P P' 

M^^, ' blocks computed in the standard formalism and the M^^, ' blocks equal to the M^^, ' blocks derived 

EB P P' 

in the pure formalism. For the remaining M,,, ' blocks, we found 



{E) (B)^ 



yrEB.BB 



i\/rEB,EB 

^^J- DDI 



9jm" -^'^ s=0 



167r 

2f + 1 
167r 



{E)i. (B) 
'^Q,l"7n"'^2-s,l"m" 



I'm" 



' + ) 7( + ) _ 7(-) 7(-) 



Re 



{E) {E)i. 



(B13) 



As is the case for the pure formalism, the EB-EE and EB-BB couplings are not zero if the spin-1 and spin-2 window 
functions have a non- vanishing ^-component. 



3. Mode- mode coupling for TB pseudo-C^ 

The temperature pseudo-multipoles are related to the CMB multipoles via the following coupling [27l |33l [34] : 

Tem = -J2Il (-1) 



(2£ + l)(2f + 1)(2£" + 1) / e e t ^ ^ ^ 



' m' i"m 



(B14) 



Combining the above pseudo-a^ with pseudo-a^^ allows us to compute the mode-mode coupling matrices for the 
pseudo-power spectra in the different formalism. The procedure is similar, though simpler, to the polarization case 
described in the previous appendix and we only list here our final results. 

First, in such the standard formalism, the mode- mode coupling matrices read 



^jTE.TE _ j^jTB,TB _ + 1 



(B15) 



Second, in the pure formalism, the mode-mode coupling matrices are given by 



njTE^TE 



lyjTB.TB 2^' + 1 



Stt 



(E) {E)i. 



'm" s 



(E) (B)^ 



01"m"^2-sl"m" 



(B16) 



Finally, in the hybrid formalism, the mode-mode coupling matrices read 



Stt 



e"m" 



0, 



(B17) 



and 



Stt 



^^J- DDI 



Stt 



(E) {E)i. 
'^Ql"7n"'^2-sl"ni" 



(E) (B)^ 
^0,£"'m"^2-s,i"'m" 



t"m" s 



(B18) 
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Appendix C: Noise bias 

The noise bias in the pseudo-C^ estimators differs from one formahsm to another. To ihustrate this point, we 
provide here the exphcit formulae for such biases assuming that noise is uncorrelated from pixel to pixel but allowing 
correlation between the different measured Stokes parameters, i.e. 

{Ns{n)Ns^{n^)) = als'{n)6{n - n^) 

where S' stands for the Stokes parameters, T, Q or [/, and assuming autocorrelation of a given map. To perform 
such a computation, we remind that [43^ 

A ^ 

^ ^ ,F/„(n,),,r,^(n,) = (-l)^-^'i?t,(a,/3,7)e-2'^^ (CI) 

m=—£ 

with a, p and 7 the three Euler angles defined by the rotation multiplication R{a, ^5, —7) = R{Oi^ (pi^ 0)R~^{Oj,(pj^O). 
For Hi = Hj, those three Euler angles are equal to zero. The ^/ stands for the Wigner rotation matrices and they 
satisfy 1)^^/(0,0,0) = Ss^s' (see [44]). With such an hypothesis about the noise properties and making use of the 
above addition theorem for spin- weighted spherical harmonics, it is straightforward to derive the noise bias in the 
three formalisms, as presented below. 

In the standard formalism first, all the noise biases vanish except for the EE and BB autospectra. For such spectra, 
the noise contribution reads: 

AffE = Mr = ^ £ (4q + ^hu) W dn, (C2) 

and 

Second, in the pure formalism, the noise biases for autospectra are given by 

Mf^ = Mr = 11^ + 4a) (w' + 4^ + ^ dn, (C4) 

while for cross-spectra they reads 

Mr = [ ^{'J^tqMW2]+4uMW2]) dn, (C5) 

~hl ^2 ^^^Q^"'!^^] - a^TuMW2]) dn. 



M 



TB 



Finally, in the hybrid formalism, the noise biases are easily deduced from the standard and pure formalism compu- 
tation to be given by 



47r 



1 /■ I 1 \ (C6) 



and 



Mr = 0, 
Mr = 0, 



= -i^l^^^ {4qMW2] - 4uMW2]) dn. 



(C7) 



Unlike the standard calculation, the pure computation induces some noise bias for the two types of odd-parity 
power spectra. However, the EB noise bias can be set to zero by switching to an hybrid approach. Moreover, the TB 
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and EB noise bias will also vanish if the noise is uncorrelated from one Stokes parameter to another or by invoking 
cross-spectra between maps coming from different experiments or different, and thus uncorrelated, detectors. 

We also stress out that all the noise biases computed above are indeed invariant under any rotation of the polarization 
reference basis. This can be straighforwardly checked by making use of the spin properties of the Stokes parameters 
and of the three spin- weighted window functions under a rotation by an angle a: 

T -^T' = T, 
(0 ± iU) {Q' ± iU') = e±2^"(Q =b iU), 
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